Irreversibility in a simple reversible model 
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This paper studies a parametrized family of familiar generalized baker maps, viewed as simple 
models of time-reversible evolution. Mapping the unit square onto itself, the maps are partly 
contracting and partly expanding, but they preserve the global measure of the definition domain. 
They possess periodic orbits of any period, and all maps of the set have attractors with well defined 
structure. The explicit construction of the attractors is described and their structure is studied in 
detail. There is a precise sense in which one can speak about absolute age of a state, regardless of 
whether the latter is applied to a single point, a set of points, or a distribution function. One can 
then view the whole trajectory as a set of past, present and future states. This viewpoint is then 
applied to show that it is impossible to define a priori states with very large "negative age" . Such 
states can be defined only a posteriori. This gives precise sense to irreversibility — or the "arrow 
of time" — in these time-reversible maps, and is suggested as an explanation of the second law of 
thermodynamics also for some realistic physical systems. 

PACS numbers: 05.70.Ln, 11.30.Er, 05.45.Pq 
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I. INTRODUCTION 

This paper is devoted to an analysis of one aspect of 
the second law of thermodynamics, namely its statement 
about irreversibility. However, I am not going to prove 
essentially new facts concerning realistic physical sys- 
tems. What I actually want to do is to illustrate the 
origin of the irreversible behavior of time-reversible sys- 
tems using an extremely simple model. The reason for 
such an approach is that, in my view, the second law (at 
least in its most general formulations) does not formu- 
late new facts about dynamical systems, it just describes 
their properties in a new form — that is why it is being 
derived from the underlying dynamics. One could say 
that the problems usually associated with substantiation 
of this law are not so much of a physical as of a concep- 
tual nature. In such a context it may be acceptable to 
rely on demonstrations instead of proofs, since the topic 
is more a question of semantics and interpretation. 

Microscopic laws of molecular dynamics are invariant 
with respect to the time reversal: only the momenta 
change their signs upon the transformation t — > — t. This 
implies that if these laws allow some evolution of a sys- 
tem, they allow also an evolution in which the system 
passes through the same spatial configurations as the 
original ones but in the reversed order (and with reversed 
velocities). If the same systems are viewed macroscopi- 
cally, they evolve, on the contrary, in one direction only: 
they demonstrate irreversibility which is formulated in 
the second law. According to the latter, only evolution 
leading to growth or preservation of disorder is observ- 
able, or in other words, there is an arrow of time. This 
conceptual asymmetry represents a fundamental prob- 



lem, called the problem of irreversibility — known for 
more than a century — which does not seem to be com- 
pletely solved up to the present time. 

There is a plethora of theories trying to explain the 
origin of macroscopic irreversibility. The approaches can 
be roughly subdivided into those treating ensembles (for 
overview see e. g. and references thereinjand the oth- 
ers studying individual systems (see, e. g. [2|)- 

In isolated systems, the irreversibility is usually being 
reduced to asymmetry in possibilities to prepare initial 
states which would evolve to equilibrium as compared to 
those evolving away from it [3|, U| . In open systems, the 
classical approach is to view the environment as a source 
of random perturbations 0, 0| , thus actually substitut- 
ing deterministic systems by stochastic ones. In such 
a way, however, the most appealing aspect of the prob- 
lem — the reconciliation of microscopic reversibility with 
macroscopic irreversibility — is being lost. 

Recently, a new promising approach to solution of the 
problem has been undertaken, studying, among other 
things, the simple model called the rotated baker map Q, 
defined on the unit square in the following way: 



B r (x,y) 



_ / (2x/3, 3y) 



for y < 1/3 



(x/3 + 2/3, 3y/2-l/2) for other y. 



This model generalizes the well-known "classical" baker 
map, 



C(x,y) = 



_ / (2s, y/2) 



(2a;-l, y/2 + 1/2) 



for0<x<l/2 (2) 
for other x, 
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and has many interesting properties. 

In this context, we will be interested in further general- 
ization of Eq. Jp) , which represents a very simple model 
enabling to demonstrate many features of interrelation 
between reversibility and irreversibility. 
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In the following discussion, we will be frequently en- 
countering the notions of measure and dimension, so that 
I find it meaningful here to give a brief description of 
what is meant by them. 

Avoiding the use of excessively technical language, we 
can say that measure /i on R™ is a set function assigning 
a non-negative number to any subset of R™ in such a way 
that the measure of the empty set is zero, the measure of 
a subset A C B is ^(^4) < n{B), and the measure of the 
union of subsets is less than or equal to the sum of the 
measures (the strict equality holding only in the case of 
disjoint subsets). This notion is a generalization of that 
of area or volume and it frequently reduces to it. 

A specific kind of measure is the dimension of a set, 
which is again a generalization of our intuitive notion 
concerning physical (topological) dimension. There are 
many ways to define and calculate it, leading to different 
values. The standard set of dimensions is based on 
partitioning the phase space into a finite number N(e) 
of disjoint e-cells (boxes) and considering the probability 
Pi (e) (the so-called natural measure) of finding points of 
the set in the box i. Of course, the probability has to be 
normalized so that 



AT(e) 

X>(«o = i. 

The dimension is then calculated, for any integer s > 0, 
according to the formula 



D x 



J-liml^M. 
s - 1 e^o In e 



(3) 



where 



N(e) 

I(s,e)=J2Pi 



i=i 



Specifically, for s — ► we obtain the Hausdorff (or 
box-counting) dimension 



Dq = — lim 



ln[JV(e)] 



+o hie 

and for s —* 1 the information dimension 



(4) 



lne 



Dx = lim 

e->0 



(5) 



which is commonly used to describe basic properties of 
fractals. In effect, the information dimension of a set of 
points gives crude information about its inhomogcncity 
(hence the name). 

It may be worth mentioning that for s — > 2 we ob- 
tain the correlation dimension. One can also prove under 
rather general conditions that D q < D r for q < r. 



II. DEFINITIONS OF REVERSIBILITY 

To introduce shortly the notion of (ir)reversibility, con- 
sider a continuous dynamical system (a flow) defined in 
a phase space T and described by differential equations. 
If St denotes the evolution operator, taking the present 
state 70 € r to the future one, 74 = St 70, then the present 
state can be retraced into the past as well. If this can be 
achieved by applying a well-defined operator S-t (origi- 
nating in differential equations), i. e., S-t"fo — 7— 1> then 
reversibility of the dynamics means that S-t is defined. 
However, the inverse operator S-t may differ significantly 
from St, so that the dynamics of the inverse evolution 
may be different from that of the forward one. If, how- 
ever, S-t and St differ just in the sign of the parameter 
t, such reversibility reduces in a functional analytic ap- 
proach to the statement that the family of operators {St} 
represents a group, defined by a generator, the latter be- 
ing closely related to the differential equations describing 
the dynamics 0. 

Another definition is not so much concerned with 
changing the direction of evolution as it is with the pos- 
sibility to "reverse" the final state 74 to a reversed one 
Tjt (typically by reversing the directions of all momenta) 
and with application of St to the latter. The dynamics 
is then said to be time-reversible if 



S t T lt = Tjq for all t > 0. 



(6) 



The question about the existence of S-t is not so impor- 
tant here and therefore such an approach seems to be 
more general than the previous one. 

Obviously, the properties of the reversal operator T 
will depend on St in general. We expect to have different 
operators T for different evolution operators St- There- 
fore, T can be any transformation which transforms a 
final state of evolution into the initial state for the same 
evolution — fulfilling, of course, the property JBJl. Every 
such transformation will be evidently idempotent, i. e., 
T 2 = I, and therefore T _1 = T. 



1-k 
k 



T<y- 



7o 
i 



7n 



T-y n 



FIG. 
ity. 



1: Commutativity as the defining property of reversibil- 



In real physical systems, the operator T has the evident 
meaning of the change of momenta: pi — > —pi. The fact 
that for time-reversible systems we have StTjt — T~/q, 
is just the consequence of the physical laws there. If, 
however, we want to speak about reversibility of systems 
in which there is no analogy to momenta — as is, e. g., 
the case of two-dimensional maps which will interest us 
in the following — we have to accept this consequence 
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as a definition and denote as reversible all maps M, for 
which there exists an operator T, making the diagram in 
Fig. H commutative, i. e., such that M n Tj n = T70, or 
more generally 

M k T ln = T 7 „_ fc . (7) 

The operator T is defined here not by the physical 
essence of M, but just by the latter requirement, and we 
have to find it. If it does not exist, the map M is not 
reversible. For reversible M, the inverse map M~ 1 will 
then be, evidently, M^ 1 = T~ X MT = TMT. 

III. GENERALIZED BAKER MAP 

We will generalize the classical map to what may be 
called "generalized" baker map B w (GBM for short), in a 
way similar to the one described in || and denoted there 
as the "slightly generalized baker's transformation" . The 
map is defined |25j for any w > 1 and is acting on points 
7 = (x, y) of the unit square E = [0, 1] x [0, 1]. For the 
points < x < (w — its action is: 

B w {x,y)=(-^-x^\ (8) 
\w — 1 W J 

and for the remaining ones (w — l)/w < x < 1: 

B w {x,y) = L(i-1) + 1, — y + i ) . (9) 
\ w w J 

Later it will be advantageous to have the above defini- 
tions rewritten in the form 

r> /_ ,a _ / for < x < (w - l)/w 

a w \x,y) - ^ ( RxX ^ Ryy } for ( w „ i)/^ < a; < l ; 

(10) 

with the evident expressions for the "left" and "right" 
operators: L x x — wx/(w — 1), R x x — wx — w + 1, L y y 
— y/w and = (wy — y + l)/w. 




w-t 



FIG. 2: B u , transforms the filled rectangle with the area 
(w — 1) /w into the rectangle with the area 1/w. Similarly, the 
complementary empty rectangle with the area 1/w is trans- 
formed into the one with the area (w — l)/w. 

B w is a piecewise linear mapping (see Fig. [SJ behav- 
ing differently to the left and to the right of the verti- 
cal line x = (w — l)/w, which we will call the dividing 



line. One can easily check that in order for the map be 
time-reversible, i. e., fulfill the time reversal condition 
B^T^in = T70, one has to define the reversal operator T 
as "rotation" around (or reflection with respect to) the 
second diagonal y = 1 

T(x,y) = (l-y,l-x) (11) 

The map can be analyzed using modern computer al- 
gebra systems, which enable us to compute the action of 
B w with absolute precision. One can then observe time- 
reversed evolution on a computer screen, and analyze the 
reversibility. However, to prevent rounding errors, such 
computations require using rational coordinates and ra- 
tional w, instead of finite-precision decimal values 10]. 

All simulations described in the paper were performed 
under the above conditions. The analysis was further 
simplified by restricting w to integers. The latter has no 
effect on the results presented in the paper, so that in 
the following I will mostly suppose integer values of w. 

The expansions caused by B w in the x direction and 
contractions in the y direction are characterized by local 
logarithmic rates l x and l y 

dB w (x,y) dB w (x,y) 

lx = In 5 and l y = In . (12) 

dx dy 

Since the action of B w is different for points lying to the 
left and to the right of the dividing line, we will have two 
rates in the x direction 

£ = In— — and l?(w)=\nw (13) 

x iu — 1 x 

and two in the y direction 

^ = In - and l*(w) = In ^— ^. (14) 
w w 

The map B w represents probably the simplest possible 
model exhibiting "microscopic" reversibility and "macro- 
scopic" irreversibility, which strongly motivates its study 
as that of an example illustrating the foundations of irre- 
versible thermodynamics, see, e. g., 11] . To start with, 
we mention that all the points of E can be subdivided 
into fixed points, cycles (periodic orbits), and attractors, 
as well as points approaching those sets. Let us consider 
each of the sets separately. 

A. Fixed points 

The map B w possesses two hyperbolic fixed points, 
(0,0) and (1, 1). The local stable manifold Wf oc for 
the point (1,1) contains any subset of the vertical line x = 
1 in E, since B™^ approaches (1, 1) for 70 = (1, yo) with 
any yo within this subset. Similarly, the local unstable 
manifold for the point (0, 0) contains any subset of the 
horizontal line y = in E (which we shall denote as 
the primary line in view of its later role) contained in 
(0, w — 1/w), because B™ 70 departs from (0, 0) for 70 = 
(xq,0) with any xq within this subset. We will discuss 
corresponding global manifolds W s and W u later. 



4 



B. Cycles 

The generalized baker map appears to have a rich 
structure of periodic orbits, or cycles. Consider a point 
7 = [x, y). Any combination of operators L x , R x acting 
consecutively on x will yield an expression linear in i, so 
that setting it equal to x will give us an equation with a 
unique solution. As an example, the equation L x R x L x x 
— x leads to 

w 3 x — w(w — l) 2 

(w - l) 2 =Xl 

with the unique solution 

w(w - l) 2 



X\ 



w 3 — (w — l) 2 



(15) 



We can find the solution for the second coordinate y\ 
similarly, using the equation L y R y L y y = y, 



yi 



w 3 — w + 1 



(16) 



The point 71 = (xi,yi) is then one of the points of the 
three-cycle (remaining two points can be calculated by 
two applications of B w ). Remembering that the opera- 
tor L x (Rx) acts on the point to the left (right) of the 
dividing line, we see that the prescribed succession of op- 
erators (we can call it the operator structure of a cycle) 
will create a cycle which will visit corresponding sides of 
the dividing line. 

There are only two structures which do not generate 
cycles, namely L n and R n , leading to stationary points 
(0,0) or (1, 1), respectively. But for any other structure 
there exists precisely one cycle visiting both sides of the 
dividing line in the given order. This shows that there 
is a one-to-one correspondence between the x coordinate 
of arbitrary point on the cycle and the sequence of po- 
sitions of remaining points with respect to the dividing 
line. Therefore, the x coordinate of a point of the cycle 
"encodes" the sequence of operators L x and R x (opera- 
tor structure), and vice versa, in a way very similar to 
that of the Bernoulli shift This is related to what is 
usually referred to as "symbolic dynamics" ^l| . 

The set of all possible cycles has interesting structure 
and deserves a detailed study. Here I just mention some 
relevant results, starting with the remark that due to 
contractivity of B w in the y direction, every cycle is a 
limit set for a subset of points in E. 

The number of cycles grows exponentially with the cy- 
cle length. There are 2 P possible combinations of two 
operators L x and having the length p. Two of them 
correspond to fixed points, and in the case in which p is 
not prime, some of the combinations (denote their num- 
ber r) may be further reduced to ones corresponding to 
shorter periods. Every p of the remaining 2 P — r— 2 combi- 
nations represents cyclic permutations, so that they cor- 
respond to the same cycle. The total number of different 
cycles of the length p is then (2 P — r — 2) /p. 



For any value of w we have, consequently, an infinite 
number of all possible cycles. The cycles created by the 
same sequences of L x and R x are topologically indepen- 
dent of w, and with growing w they are just scaled to 
smaller dimensions (the denominators are growing) and 
translated towards the point (1, 0). This implies that 
there is the same number of cycles for any w. In the 
case of w = 2, every point with an odd denominator 
of the x coordinate belongs to some cycle. For w > 2, 
however, we do not have such a simple rule due to the 
above-mentioned scaling. 

Consider now a period with very large p. Only a very 
small proportion out of 2 P sequences of corresponding 
operators L x and R x will be ordered in any sense - 
the majority will look like random sequences. The orbits 
they will generate will therefore be indistinguishable from 
random sets of points on the trajectories. 

It is evident that the set of points belonging to cycles, 
and of points approaching them, is of zero measure in E 
and consequently the behavior of cycles is not what we 
could observe in GBM frequently. 

We now come to the most important subset of E, which 
is — as we shall see — dense in E and therefore the be- 
havior of its points represents typical properties of GBM. 



C. Attractors and their properties 

Averaging the logarithmic expansion rates i|13|) 
and l|14f> over typical trajectories in E (see ^3$), one •ob- 
tains what is usually called Lyapunov, or time-averaged, 
exponents — the positive one 



W — 1 w 
Ai = In 



1 lnu; = — In 

w — 1 w w 



(w - l) w - 1 (1 ' 



and the negative one 



. w — 1 1 1 . ui-1 1 . tii-1 ,,„. 
A 2 = ln — + — ln = — ln . (18) 

w w w w w w w 

From the above expressions one sees that with growing 
w > 2, both exponents are monotonously decreasing, 
and their limit behavior for w — > 00 is X\(w) — > and 
^{w) — > -00. 

The existence of the positive Lyapunov exponent sug- 
gests that one should expect chaotic behavior of the it- 
erates of -B^7o (for almost every 70 S E), and the ex- 
istence of the negative one the existence of a (strange) 
attractor Both are observed when one iterates B w 
beginning with almost any starting point — see Fig. |21 

It follows from previous remarks that exceptions to this 
general statement include points of Wf oc , approaching 
the hyperbolic fixed point (1, 1), and points belonging to 
cycles, or approaching them. 

One can prove that the attractor consists of an infinite 
set of lines parallel to the x coordinate, see Fig. 0] and 
can be generated by successive applications of B w to the 
primary line. The construction is based on the iterated 



5 



function system, see [H| . It is well known that in the case 
of Xi + A 2 < 0, the attractor is a (multi)fractal object. 
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FIG. 3: Iteration of the point (95/100, 95/100) by the action 
of Bs discloses a distinct attractor with clearly visible self- 
similarity. 10 000 iterations are shown. 



From the definition of an unstable manifold, it follows 
that the points lying exactly on the attractor (i. e., points 
generated by the described construction) represent the 
global unstable manifold W u for the point (0,0). In this 
sense, the latter represents the source for the attractor. 

The attractor is markedly inhomogeneous (strange) for 
greater values of w. With growing y, its density decreases 
(for w > 2, but increases for w < 2 if we permit rational 
non-integer values for w). With decreasing w > 2, the 
inhomogeneity is less and less pronounced, until at last, 
for w — 2, the limit object becomes the set of equidistant 
horizontal lines. The latter specific case will be treated 
separately in Sec. IIXI 



w 


2 


3 


4 


5 


10 


100 


1000 


Di(w) 


2.000 


1.734 


1.506 


1.376 


1.156 


1.012 


1.001 



TABLE I: Information dimension of the attractor of B n 
some values of parameter w. 



for 



As with any attractor, one would like to know its di- 
mensions £) s (w), s > 0. The Hausdorff dimension is 
evidently Dq(w) = 2 (see e. g. 14]). The information 
dimension is calculated with the help of formula (J5J , tak- 
ing into account that one can use, in place of probability, 
the density of iterates of a typical single point yal . Re- 
peating the calculations of Hoover and Posch [13j. but 
for general w, one arrives at the following expression for 
the information dimension J5j|: 



I>iH = 1 



lnK; {w - iy- w ) 

ln((io - l)w- w ) 



(19) 




4 5 6 7 
generation 

FIG. 4: Iterative generation of the attractor of B w (w = 5). 
First, B w is applied to the line segment y = (primary line). 
We call the resulting two lines (including the primary line) 
the first generation. Applying B w repeatedly to all lines of 
the previous generation, we obtain 2 n prefractal lines after 
n iterations. New generations are illustrated by gradually 
shorter lines. 



Table I gives the dimension for a few values of w, and 
the graph in Fig. [S] enables us to see the overall depen- 
dence of D\(w) on w. Evidently, one can obtain any 
value of Di(w) from the interval 2 > Di(w) > 1, by con- 
trolling uu. Since D\ can be viewed as the measure of 
attractor inhomogeneity, the approach of the dimension 
to the value of 1, with growing w, suggests that the it- 
erated points tend to accumulate (condense) on smaller 
subsets of E. 




FIG. 5: Functional dependence of information dimension 
Di(w) on the value of parameter w. One sees that the dimen- 
sion will approach 1 with infinite growth of w. The behavior 
for 1 < w < 2 does not represent anything new since for these 
values the unit square E is just contracted in the direction of 
growing y. 



It can be proved [l5j that every attractor generated by 
an iterated function system is the closure of its periodic 
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points. This suggests that the evolution generated by B. 
may show Poincare recurrences. 



IV. REPELLORS AND SYMMETRY OF 
EVOLUTION 

Since our main preoccupation is the study of interre- 
lation between reversibility and irreversibility we cannot 
be interested only in the "future" of a state defined by 
an initial point 70, i. e., in the trajectory beginning in 
this point, but we have to consider its "past" as well, 
viz. the trajectory ending there. Here we will therefore 
try to find out what is possible to tell about the whole 
trajectory going through a given point 70. 

Repeated application of B w to arbitrary point 70 gen- 
erates a sequence of future points {71, 72, . . . 7n-i, 7n}, 
where 7^ = -B£,7o- Let us denote by j-k a point from 
which the point 70 ensued after the application of the 
map B^. Then the past of the point 70 will be repre- 
sented by the sequence {7_„,7_„ + i, . . .7-2, 7-i}- Evi- 
dently, 7_ fc = TB k T~f , and we will call the inverted 
sequence {7-1, 7-2, ■ ■ • , 7-n+i,7-n} the backward itera- 
tion of the state 70. The fact that it can be obtained also 
by repeated application of B^ 1 to 70, is irrelevant here. 

It is clear that T70 is arbitrary (or random) in exactly 
the same sense as 70, so that the iterates B^T^o will 
approach the attractor as well. This has the following 
consequences. Since B w is contracting in the y direction, 
the distance between any two points, having exactly the 
same x coordinate, will quickly decrease under the action 
of , and the properties of the future trajectory will be 
determined essentially by the x coordinate of the initial 
point. The reversal T interchanges the components of 
coordinates, so that during the backward iteration the y 
coordinate of the original point will similarly determine 
the past trajectory. We can thus say that the information 
about the global aspects of the future is contained in the 
x coordinate of 70 and the information about the past 
in its y coordinate. We will express this fact by writing 
symbolically 7 = (x futl y past ). 

This, however, means that the past of a point will ap- 
proach (under backward iteration) an object symmetrical 
with respect to the attractor. As this object is composed 
of vertical lines, every point in its vicinity will move away 
from it under forward iteration, since the differences in 
x coordinates are growing with iterations. That is why 
it is being called the repellor. 

We have thus come to an important result: an or- 
bit going through arbitrarily chosen point 70 has its past 
close to the repellor, and its future close to the attrac- 
tor. In other words, a typical trajectory generated by 
the action of B w departs (in the past) from the repellor 
and approaches (in the future) the attractor. Such unidi- 
rectional behavior is sometimes characterized as demon- 
strating the existence of the "arrow of time " and here it 
is the direct consequence of the dynamics defined by B w . 
The time reversal of such a trajectory will again be go- 



ing from the repellor to the attractor, since T transforms 
future points into the past ones (and attractors into re- 
pellors) and vice versa. The time reversal is thus not able 
to change this global aspect of evolution. 

1-1 1 




0.2 0.4 x 0.6 0.8 1 



FIG. 6: Symmetry between the past and the future of the tra- 
jectory (for w — 5/2) with the initial point So = (1/10, 9/10) 

— denoted by a box — lying exactly on the diagonal. 

This state of affairs can best be seen when we choose 
a point So lying exactly on the second diagonal, i. e., a 
point Sq = (xq,1 — Xq). Such an initial point remains 
unchanged upon reversal and therefore its past will be 
the reversal of its future: past and future parts of the 
trajectory will be exactly symmetric with respect to the 
diagonal, S-k = TSk, see Fig. Or still in other words, 
we obtain the past part of the trajectory, in this specific 
case, by reversing its future. 

Considering the explicit construction of the attractor, 
described in the previous section, we immediately see 
that the points lying exactly on the repellor represent 
the global stable manifold W s for the point (1,1). In 
this sense, the latter represents the sink of the repellor. 

Returning to cycles, we discover an interesting prop- 
erty, applying to any w. There are trajectories which un- 
wind from a p\ cycle in the past and approach a pi cycle 
in the future; see the example in Fig.0 In my view, such 
behavior illustrates best the discussed difference between 
the past and the future, since it demonstrates that there 
exist two different limits in the behavior of B w , separated 
by intermediate states. Let us mention in passing that 
the period of the future and the past cycle is determined 

— in accordance with the general rule expressed above 
symbolically as 7 = (xf u t,y pas t) — by denominators of 
the coordinates x and y, respectively. 

B w can be viewed as mimicking the behavior of a ther- 
modynamic system (a nonequilibrium thermostatted one 
for w 7^ 2, and an isolated one for w = 2) and therefore 
it could be employed to explain the origin of irreversibil- 
ity. In this way, one is lead to the tentative conclusion 
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u 0.2 0.4 x 0.6 0.8 1 

FIG. 7: Trajectory unwinding from a 6-cycle in the past and 
approaching a 4-cycle in the future (for w = 3), "generated" 
by the initial point (16/73, 63/103), denoted by a box. 



that the evolution of some physically relevant systems, 
sharing common properties with B w , is such that their 
states approach an attractor and in the past they de- 
part from a repellor. This being the case, we arrive at 
an alternative explanation of the interrelation between 
reversibility and irreversibility, avoiding some paradoxes 
usually derived from the second law. For example, the 
question of why we do not observe the evolution going 
in the direction from attractor to repellor now becomes 
meaningless, because such a possibility is ruled out by 
the system's dynamics. The problem is instead being 
shifted to the question of why we encounter an approach 
to the attractor but not the departure from the repel- 
lor — both proceeding in one and the same direction of 
time. In other words, the question is, why do we observe 
only one part of the full (and in principle possible) evo- 
lutionary trajectory? The next sections are devoted to 
tackling the problem of irreversibility in this alternative 
formulation. 



V. AGE OF STATES 

Since the typical full phase trajectory (considered un- 
bounded both in the future and in the past) consists of 
points moving from the past repellor to the future attrac- 
tor, it is quite natural to call the points between the past 
and the future the present ones. This leads us to con- 
sider the possibility of introducing an "age" that could 
be applied along the trajectory. 

Looking at the definition (JSJ and © of B w with ratio- 
nal w, we see that if we choose arbitrary point 70 with 
rational coordinates, the coordinates of 7„ = B™jq will 
be rational too, and with growing n they will be repre- 



sented as fractions of growing integers (except the cases 
of periodic orbits) — they will become more "complex" . 
The same applies to backward iteration 7_fc = TB k T^. 
Since, by definition, the denominator of any coordinate 
in E is not less than its numerator, we can assess the 
"complexity" or "simplicity" of coordinates by the value 
of denominators: smaller denominators will mean "sim- 
pler" coordinates. 

Expressing an arbitrary point on the trajectory (us- 
ing fractional expressions for rational coordinates) as 
7 = (x,y) = (p/q,r/s), we realize that in the (distant) 
future the mapping of the y coordinates of points is dom- 
inated by y — > y/w, whereas that of the x coordinates is 
dominated by x — > x/(w — 1). This means that the y 
coordinates will have greater denominators than the x 
coordinates (s > q), whereas in the (distant) past it will 
be the opposite, since the role of coordinates becomes ex- 
changed. We could thus base the notion of "age" on the 
difference s — q. For future points the latter will be posi- 
tive; for the past ones, negative. However, except for the 
points lying exactly on the diagonal (see Fig.|fJJl, the pres- 
ence will not be sharp: we will not observe monotonous 
behavior in its neighborhood and different trajectories 
may have present states of different ages. This is espe- 
cially true if we choose present points with q^> s, making 
them seem to be past in such a way. 

To illustrate the situation, consider the following typi- 
cal subsequence of values s — q for w = 3 and k = — 8 . . . 8, 
with initial point 70 = (1/4,3/5): 

. . . , -2027, -649, -203, -61, -17, -8, -2, 1, 1, -3, -1, 13, 
103, 341, 1151, 3517, 10679, .... 

The subsequence is centered around the value 1, denoted 
in bold, corresponding to the initial point 70. The differ- 
ence s — q lacks the most important aspect of age, namely 
its linearity. We can obtain the latter by defining the age 

as 

r = sgrx{s-q)\og w {\s-q\) (20) 

and skipping eventual (and rare) zeros for which this ex- 
pression is not defined. In such a way, we get a really very 
precise approximation to linear age — a result which can 
be understood by noticing that B w always maps s — + ws, 
whereas q either remains unchanged or q — ► (w — l)q. 
Then for 7fc we have on average s — q « w k — (w — l) fc_1 
for k > 0. The log w of the latter approaches, with grow- 
ing fc, the value of k+const with very high precision. The 
above sample (sub)sequence is then transformed into the 
following one: 

. . . , -6.93, -5.89, -4.84, -3.74, -2.58, -1.89, -.63, 0, 0, -1, 
0, 2.33, 4.21, 5.30, 6.41, 7.43, 8.44, . . . 

which approaches, with growing k, equidistant values of 
"age" , separated by exactly unit steps. One has to note 
that corresponding points on different trajectories can 
have slightly different ages. 

These reasonings demonstrate that in a microscopi- 
cally described system of GBM, the age of any point can 
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be found in an unambiguous way and with reasonable 
precision; see Fig El It will increase in equal steps from 
past to future — except the small neighborhood of the 
"present" state. This age is absolute in the sense that it 
does depend only on the distance from the present, i. e., 
on the number of iteration steps. This is due to the fact 
that coordinates of points contain "traces" of their evo- 
lution, which enable us, at least in principle, to recover 
the iteration index k from the coordinates of arbitrary 
point. It is appropriate to stress here that the case of 
the classical baker map (w = 2) is no exception to this 
general rule. 
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FIG. 8: Dependence of "absolute age" r on iteration in- 
dex k. Although one of the worst cases, with initial point 
(50 001/100 000, 11/30), has been chosen here intentionally, 
the age r differs from k by less than 10. This situation is 
far from typical, because in a random choice of a point it is 
not probable to obtain such a great difference in the preci- 
sion of its two coordinates. The typical behavior is character- 
ized, contrary to the depicted one, by two parallels close to 
r(fc) = k, exhibiting a sudden jump around k = 0. 



The described approach to age might seem extremely 
artificial, but it can be viewed, nevertheless, as represent- 
ing what may be called "microscopic" age. Of course, one 
should not overemphasize this model description, but nei- 
ther should it be underestimated. In my view, it demon- 
strates — together with the above-mentioned behavior of 
cycles, illustrated in Fig.0 — that there is, at least qual- 
itatively, a difference between past and future already at 
the microscopic level. This then implies the existence of 
the property of the system, which I would call micro- 
scopic irreversibility. 

If we now consider the evolution of a set of N points 
under the action of B w , the above reasonings will apply 
to any point of the set. We could therefore — at least 
in principle — define in any evolutionary trajectory of 
the set the "age" as the average of ages of all individual 
points. The "present" state would then be the one with 
the (averaged) age closest to zero. 



However, looking at the set without the knowledge of 
the coordinates of individual points, it is possible to know 
its "age" also by observing its shape. Specifically, in the 
case of a subset of E having a simple boundary (e. g., 
the case of points lying within a small square), the B w — 
causing the growth of fractional expression of coordinates 
of all the points — will cause also the "distortion" of 
the boundary. This then means that such distortion is 
commensurate with the age of the set. To make this idea 
more definite, we turn to the distribution functions on 
E, giving the density (of probability) of points in E. 




FIG. 9: Two states in the evolution of a distribution function 
(under the mapping -B3), the present state of which is defined 
as constant on the unit square. 



The map defined by B w on E will translate in a 
straightforward way into the evolution of distribution 
functions having E as their support. The distribution 
function, constant on a compact subset of E, will in the 
forward evolution approach a function independent of x, 
and in the backward one, independent of y (see Fig. |5J|. 
But even general distribution functions will very quickly 
approach functions almost independent of a;, or y. That 
is why one can characterize, with sufficient precision, the 
future limit behavior of any distribution function with 
the help of the function f(xo,y), representing a section 
across f{x,y) by a plane perpendicular to the x axis, go- 
ing through arbitrary point xq, see Fig. 1101 For such a 
function of one variable, a variation on the unit interval 
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is defined as the supremum 

m 

V y (f) = S upJ2 \f(xo,Vj) - f(xo, Vj-i)\, (21) 

3=1 

taken over all possible dissections of the interval by m > 1 
points. We can regard this expression as an acceptable 
approximation of the variation of the distribution func- 
tion along the y axis. 
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FIG. 10: Section across the distribution function 0s = B| <f>o, 
with the initial function 4>o(x,y) = 1, by a plane perpendicu- 
lar to the x axis. The maximum value is 4>5 (0, y) = 32. 



In quite a similar way, one can characterize the past 
limit behavior of any distribution function with the help 
of the variation V x (f) of the function f(x,yo), represent- 
ing a section going through arbitrary point y . The dif- 
ference V y (f) — V x (f) of variations along y and x will 
then grow steadily for n from — oo to +oo, and we can 
use it to estimate the "age" of the state, similarly as we 
used the difference of denominators in the case of a single 
point. 

This then gives us the general picture of the evolution 
of a distribution function under the action of the map B w . 
The absolute value of its variation will grow indefinitely 
for 1 77, | — > oo, attaining somewhere in between its mini- 
mum. Choosing one of the states having low variation 
for the "present" state and denoting the corresponding 
iteration index by k — 0, we will again introduce what 
may be called the "age" k of the system in question. 
We could also go a bit further and make this age linear 
by taking the logarithm of the difference of variations, 
similarly as in the case of a single point. Such an age 
would then be mathematically well defined everywhere, 
with the possible exception of the neighborhood of the 
"initial" distribution function. 

Let us note in passing that the distribution function 
that evolved from standard "testing" function eV)(x,y) = 
1 will have the future maximum values (for w > 2) close 
to the line y = and the minimum ones close to the line 
y = 1. After n iterations, the values will be (w — 1)™ 
and (w — l) _n , respectively. This can be readily seen by 



taking into account that the area (w - to the left 
of the dividing line is contracted by the action of B w to 
the area 1/w, i. e., by the factor of w — 1. In a similar 
way, one will understand the behavior of the minimum. 

The ratio of the maximum to the minimum value is 
then (w ~ l) 2r \ which becomes — for w > 2 — a huge 
number after just a few iterations. This illustrates very 
well the growing inhomogeneity of the distribution func- 
tions under the action of GBM. 



VI. THE PROBLEM OF IRREVERSIBILITY 

Having introduced the notion of absolute age in a 
way formal ly d ifferent from the one used by the Brus- 
sels school [3 El] ! but conceptually — as it seems - 
rather close to it, we are ready to explain the appar- 
ent contradiction between microscopic reversibility and 
macroscopic irreversibility of generalized baker maps, or 
in other words, to explain the essence of the second law 
of thermodynamics as applied to this simple model. 

The so called "problem of irreversibility" is usually ex- 
pressed in two distinct formulations: one concerns the 
evolution of a state prepared by an experimenter, the 
other applies to the observation of a system without any 
question about its origin, i. e., of a system with an un- 
known past. We are thus facing two questions (both 
equivalent to the Loschmidt objection) which are being 
posed since the first explanation of the second law by 
Boltzmann in 1872: if both directions of evolution are 
microscopically equally possible, then (a) why we are not 
able to prepare states which would evolve in a direction 
"prohibited" by this law, and (b) why do we never ob- 
serve such evolution in systems with an unknown past? 
In the language of the absolute age, these questions re- 
duce to the following ones: why are we not able to pre- 
pare, nor to observe, states which belong to the distant 
past on the evolutionary trajectory of the system under 
consideration? 

We will treat the two questions separately because the 
explanations will turn out to be different for the two 
cases. We will also treat separately two different ap- 
proaches, corresponding to microscopic and macroscopic 
descriptions. 



VII. MICROSCOPIC IRREVERSIBILITY 
A. Irreversibility in evolution of a single point 

We have shown that if we choose a point (with rational 
coordinates) in E in any way, the coordinates of its future 
and past iterations will be almost always (i. e., with the 
exception of zero measure of points) more complex than 
the present ones. The claim that we can prepare only 
absolutely "present" points will then become actually a 
self-evident tautology, stating that we can prepare "sim- 
ple" coordinates more easily than the "complex" ones. 
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This tautology, however, expresses exactly the essence of 
the problem. If we would ask why it is not possible to 
prepare an "old" point, the coordinates of which would 
become simpler in the future, we could first of all point 
out that if that were even possible in principle, such sim- 
plification would last only for a finite number of iterations 
— until the coordinates would become relatively simple 
(the state would become "present"). However, to prepare 
ab initio a point close to the repellor, i. e., a point with 
specific and thus very complicated coordinates, means to 
prepare an exceptional point, and this is essentially more 
difficult than to generate a point at random. I propose to 
call this state of affairs, in the given context, microscopic 
irreversibility. 

It is evidently possible to prepare such a point, even 
numerically, only by choosing a point 70, letting it evolve 
to the point B^jo, and at last inverting the latter. The 
point TB^jo, having the age of —k iterations, would then 
approach, under the action of B w , the point T70 and the 
coordinates of its k iterations would be simplified. But it 
is not possible to find a general construction nor formula 
defining a set of such "old" points. 

From this point of view, the microscopic irreversibil- 
ity is related to substantially different requirements for 
preparation of the present and past states, respectively. 
Evolution leaves "traces" on the coordinates of the 
evolved point, and inverse evolution would require us, 
even in the case of a single point, to reconstruct precisely 
its whole history. But nothing similar is being required 
to prepare the present state. This crucial asymmetry 
lies, in my opinion, at the root of the explanation of irre- 
versibility 123. One should not, however, forget that this 
situation is the result of chaotic properties of the evolu- 
tion induced by the map B w . If the map did not lead to 
chaotic dynamics, there would be no tautology; nothing 
would be self-evident in this sense. 

We assumed that phase points have rational coordi- 
nates. The considerations of the last two paragraphs 
hold, however, equally for irrational coordinates. In par- 
ticular, the fact remains that inversion of evolution would 
require us to recover exactly its whole history. In this 
sense, the argument applies to any coordinates. 

Here we come to a variance with traditional view of ir- 
reversibility as formulated, e. g., by Bricmont [l8|: "All 
the familiar examples of irreversible behavior involve sys- 
tems with a large number of particles (degrees of free- 
dom). If one were to make a movie of the motion of 
one molecule, the backward movie would look completely 
natural." Certainly, playing back the movie of a single 
point moving under the action of GBM, we would not be 
able to distinguish the backward movie from the forward 
one. But why should the physical properties depend on 
our common way of observation or on our observational 
possibilities? And what finite number of points is suffi- 
cient to distinguish between forward and backward mo- 
tion? Superimposing all the phase points of the trajec- 
tory of a single particle onto one frame, as in Fig. [31 we 
would be able to make the distinction between the past 



and the future, between impossible and natural. 

The simple model system of N noninteracting parti- 
cles in a rectangular box 0] shows that large number 
of particles alone is not sufficient for irreversibility. The 
arguments of this section support this line of thought 
in that the dynamics of some systems may bring about 
the unidirectionality of microscopic evolution, manifest- 
ing itself in the appearance of attractors (and repellors). 
Observing then systems of N particles in the usual way, 
the inherent microscopic irreversibility becomes more and 
more pronounced with growing N, but it is present in the 
system all the time. 

Moreover, the approach connecting irreversibility with 
a large number of particles is not straightforward either: 
it requires, e. g., thermodynamic limit N — * 00 to obtain 
irreversibility. For further details, see, e. g., [2(j, giving 
a nice treatment of the difference between chaotic and 
irreversible behavior of Hamiltonian systems. 

It might be appropriate to mention at this point the 
Poincare recurrence theorem, stating that any isolated, 
locally measure-preserving system will eventually return 
close to any of its previous states. Approach to an attrac- 
tor, an object dense in E, has such a recurrent property. 
But the property of being close in phase space does not 
mean being of a similar age: in an arbitrary close neigh- 
borhood of any point, there can be points of any age. 
The mixing property of B w implies that the distance in 
the sense of metrics does not imply the distance in time, 
nor vice versa. In this sense the recurrence does not con- 
tradict irreversibility. 



B. Irreversibility in evolution of many points 

The analysis of behavior of a single point may be of 
interest on its own, but the essential properties of the 
model are seen much better when observing the action of 
GBM on a set of points. 

The initial state can be, in the case of GBM, prepared 
essentially by defining a subset of E and filling it uni- 
formly with random points (or according to an appropri- 
ate distribution function). The evolution of such a set 
can be observed in the author's animation [2^, where 
the subset can be chosen as a rectangle. Watching the 
evolution, one can see that points contained originally in 
the bounded subset will in the future get into a growing 
number of horizontal stripes, and in the past into vertical 
ones. Going from the past to the future, they will pass 
through the "present" subset. This situation is depicted 
also in Fig. 1111 

Consider a set fl^k of N points from which a set f2o, lo- 
cated within a small square, will evolve after k iterations. 
The greater the k, the closer f2_fc is, in the usual sense 
of metrics, to a set £l r of points randomly distributed 
over (the repellor in) E. However, irrespective of how 
close the latter two sets become, their future behavior 
will be essentially different. After k iterations, the first 
will become fio, the second will look like the attractor. 
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FIG. 11: Three states in evolution (with w = 4) of an initial 
set of 2500 randomly chosen phase points, concentrated in 
the square subset of E. We are readily able to distinguish 
between the past, the present, and the future. The projections 
of points on the diagonal y = 1 — x represent an analogy to 
the positions of physical particles. 



The first question pertaining to the second law can be 
formulated in this case as asking why it is not possible 
to prepare a state corresponding to a distant past. To 
answer it, it is sufficient to note that the evolution of a 
set of points, although being reversible, does not offer 
an equal possibility for preparation of all initial states. 
Present states — e. g., points filling the small square in E 
— can be prepared without the slightest problems, but to 



prepare a state n steps before would require localization 
of points in 2™ vertical stripes of E, and in every stripe 
with different, and well specified, density (cf. Fig. ITO)) . 
How impossible this is can be best understood when we 
realize that this task is equivalent to the one of preventing 
the occurrence of points in the complementary subset. 

Then the answer to the above question is that prepa- 
ration of past states would require us to create them ac- 
cording to a distribution function with wild variation, 
which is practically impossible. Moreover, this function 
is extremely close to the other one, obtained by averaging 
of the former, and describing a readily preparable present 
state. So the preparation would have to be, in addition, 
extremely precise. Here again we see the crucial role of a 
sensitive dependence on initial conditions, since two dis- 
tribution functions, differencing "macroscopically" only 
slightly, have an essentially different future. 

The second question, why we do not observe states 
typical for the distant past, can also be answered easily: 
we do not observe them because any such state quickly 
becomes "present" or future, and observation means that 
we are watching a system with an unknown past without 
influencing it. 



VIII. MACROSCOPIC IRREVERSIBILITY 

Let us find out now what we can learn about GBM if 
we want to consider it as a model of a physical system, 
exhibiting the difference between microscopic and macro- 
scopic behavior. The difference derives typically from the 
distinction between the full description in phase space 
and various reduced descriptions. In the case of GBM, 
the unit square E can be viewed as the phase space of one 
particle (denoted usually as the /i-space), and any set of 
points in it as the complete description of a state of a ID 
system of noninteracting particles. If we want to use the 
model to illustrate the micro/macro dichotomy, we have 
to find a counterpart to the less complete description. 

We depart from the inversion operator T, the effect of 
which is analogous to the inversion of momenta in classi- 
cal mechanics. If the GBM has to resemble, at least re- 
motely, the dynamics of ID gas, then the particle position 
should be represented by a variable which is an invariant 
of T. The projection of points in E onto the "second" 
diagonal, y = 1 — x, has this property and we therefore 
assume that it represents positions. The projection of the 
attractor looks like a stationary (non) equilibrium distri- 
bution of points, which can be called the "limit state. " 
Due to the symmetry with respect to the diagonal, the 
projection of the repellor will be identical to that of the 
attractor. The orthogonal projection of a point in E onto 
the "first" diagonal, y — x, could be viewed, with certain 
reservation, as the "momentum" of a particle. 

The "macroscopic" behavior of the system of particles 
is then described by what can be seen on the second di- 
agonal (see Fig. I12fl . There we no longer observe the 
approach to the attractor, nor departure from the re- 
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FIG. 12: The projection of evolution depicted in Fig. llll onto 
the second diagonal (scaled to unity length). Only 500 points 
were used to make the structure visible. The left-hand side of 
the line corresponds to the point (0, 1) in E. We are evidently 
not able to distinguish between the past and the future of the 
projected set, but we readily see the difference between both 
on the one hand, and the "present" on the other. 



pellor. In the direction from present to future, we see 
only the change of inhomogeneous distribution of posi- 
tions into the time-independent limit state. In the direc- 
tion from past to present, we see the opposite change of 
almost homogeneous (non) equilibrium state into the in- 
homogeneous present one. Moreover, we are not able to 
distinguish between the past and the future due to the 
attractor/repellor symmetry. The transition from past 
to present contradicts the second law for a gas on a line; 
the transition from present to future is in accord with it. 

We can now answer the two questions posed tradition- 
ally with respect to the second law. We already know 
that we cannot prepare microscopically the "past" state 
because this is practically impossible. But the macro- 
scopic preparation offers still fewer possibilities, since in 
this case we can "assign" only positions to particles, but 
not "momenta" . If we were even somehow able to create 
some specific arrangement of points on the diagonal, we 
would know nothing about their "momenta" , so the lat- 
ter would be necessarily random, and consequently the 
phase points would not be on the repellor. This explains 
why it is not possible to prepare macroscopically a state 
which would evolve contrary to the second law, and it 
seems to be the full explanation. 

We cannot observe evolution going from the repellor 
to the present state for different reasons. First, any ob- 
servation of a closed (isolated or thermostatted) system 
has finite duration. If we regard the trajectory generated 
by as a sequence of states with unbounded past and 
future, then the probability of observing any specific fi- 
nite sequence of states is zero. The observation of the 
present state is by definition specific and finite, so the 
probability to find it at random (on a randomly chosen 
trajectory) is vanishingly small. But the probability of 



picking at random arbitrary interval with exclusively past 
states (on the repellor) is the same as that of seeing one 
with exclusively future states (on the attractor) — i. e., 
1/2. However, we cannot distinguish macroscopically be- 
tween the past and the future states. This explains why 
we are not able to know that we have seen a system in 
the past even if it would be there. This holds even under 
the assumption that systems can exist under unchanged 
outer conditions indefinitely. 

But the trajectories of closed systems do not have un- 
bounded past. In the finitely remote past a state "came 
into existence" as a result of interaction of the system 
with the surroundings, and afterwards its boundary con- 
ditions remained unchanged. That "first" state was ac- 
tually the "present" one, and from the moment of its 
appearance the system tends to the attractor. There- 
fore, the unbounded past of the presently observed state 
is actually an illusion: from the existence of the present 
state we can assume the existence of the past states, we 
can even in principle recapture them, but the mere ex- 
istence of the present state does not guarantee the real 
existence of past states on the repellor. Thus we can find 
the system always only in the future, approaching the at- 
tractor. This explains the second law in the form stating 
the properties of observed states, because it describes the 
behavior of a system with an unknown past. 

This leads us to formulate the macroscopic irreversibil- 
ity of the GBM as meaning that the evolution of this sys- 
tem approaches an attractor, and if it would be possible 
to recover its evolution from the past, the system would 
depart from the repellor. Such formulation is not sub- 
ject to the Loschmidt paradox; it is free of any seeming 
contradictions. In my view, these are sufficient reasons 
to accept it. However, to do this, we have first to take 
into account that the usual statement of irreversibility 
applies to conservative (isolated) systems which should 
be simulated by B2 — the classical baker map. 



IX. THE CASE OF THE CLASSICAL BAKER 
MAP 

In the previous reasonings, we were dealing mostly 
with the full interval of parameter values, w > 1. How- 
ever, the specific case of w = 2, corresponding to the 
classical baker map, does not share all general proper- 
ties stated above, and the question therefore arises as to 
which of the properties do not carry over to this specific 
case and what consequences this may have. 

The attractor of the GBM reduces, for w = 2, to the set 
of horizontal lines with = Ci/2 k , where k denotes the 
generation number and Cj < 2 fc is an integer. Similarly, 
the repellor becomes the set of vertical lines with x coor- 
dinates of the same type, Xik = Ci/2 h . Speaking mathe- 
matically, these sets represent an attractor/repellor pair 
as well (see e. g., p. 82), therefore we will continue 
to denote them so. 

From the physical point of view, however, those ob- 
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jects are usually considered irrelevant, since they are not 
"strange ". All dimensions of the attractor and repel- 
lor namely reduce to 2, so that the latter are homogeneous 
and consequently do not represent readily "visible" ob- 
jects. Moreover, the sum of Lyapunov exponents 

2 — in 

Ai+A 2 = ku>-l), (22) 

w 

which should be negative for strange attractors, becomes 
now zero, which means that the measure of subsets of E is 
not being contracted under B 2 . Such a homogeneous at- 
tractor, nevertheless, "attracts" the points iterated by B 2 
in the important sense that the iterates visit the neigh- 
borhood of any point of E with equal probability. 

Returning now for a while to cycles, we have observed 
that the number of periodic orbits containing points with 
denominator not exceeding a given number d grows with 
falling w > 2, until in the case of w = 2 all points with a 
rational x coordinate, differing from those of the attrac- 
tor/repellor ones, lie on trajectories approaching cycles. 
This can be understood considering that the B x reduces 
in this specific case to x — » 2x mod 1 and the condition 
of a p cycle, B^.x — x, then has a unique solution for any 
rational x. The attractor is therefore approached only by 
points with irrational x coordinates, and as the former is 
uniformly dense in E, this behavior can be understood 
as chaotic approach to homogeneous coverage of E. But 
the approach to very long cycles is indistinguishable from 
the approach to a homogeneous attractor, so that also the 
majority of trajectories with rational points will seem to 
approach the latter. The typical observable behavior of 
all points — rational and irrational — will then look like 
chaotization, or approach to equilibrium. 

Accepting the terminology according to which the at- 
tractor (repellor) exists even in the classical baker map, 
being only "invisible" there, we can formulate the ob- 
tained results in a unified manner, treating classical baker 
map, caricaturing an isolated system approaching equi- 
librium, in the same way as thermostatted systems ap- 
proaching nonequilibrium steady states. 

The notion of absolute age has the same meaning for 
w = 2 as for w > 2. The present states are defined 
by simple coordinates, and the present distribution func- 
tions have small variation. One can tell the difference 
between past and future in the evolution of a point, and 
choosing an inhomogeneous initial (present) distribution 
function, one can distinguish microscopically between 
past and future not only by absolute age but also by 
predominantly vertical or horizontal stripes of maxima, 
respectively — at least in the initial stages of evolution, 
not too far from the present states. 

Absolute age is essential, in my view, for the expla- 
nation of irreversibility. What was told about it above 
applies equally to the specific case of w — 2 and sup- 
ports the view that what makes the difference between 
reversibility and irreversibility is actually observability: 
macroscopic irreversibility is observable, whereas micro- 
scopic is not — it is hidden for macroscopic methods of 



observation. In the case of w = 2, it is hidden also for 
microscopic methods. Only the generalization to w =/= 2 
is able to disclose its microscopic existence. And that was 
the reason for the study of the generalization of locally 
measure preserving maps. 

In this context, we should not leave unnoticed the Zer- 
melo paradox, relying on the Poincare recurrence theo- 
rem. According to it, even if a nonequilibrium state of 
an isolated system evolves to equilibrium, it must even- 
tually return close to the initial state, thus violating the 
second law. This means, as applied to B 2 , that a sys- 
tem of points, starting e. g., in a small subset of E, will 
eventually return close to it. Apart from the standard 
reply that the recurrence time is enormously long, the 
model of B 2 enables us to give further refutation of the 
objection. Since all rational points approach cycles, such 
recurrent behavior, albeit not excluded, would be possi- 
ble only for extremely exceptional sets. Choosing namely 
the same denominator of x coordinates of all points, one 
obtains — as a rule — the sets "condensing" on a limited 
number of points in E. However, it is not too difficult to 
find also points, which will approach periodically mov- 
ing sets. This confirms the existence of recurrences, but 
also demonstrates how rare and physically irrelevant they 
are [2j|. In this sense, the recurrence does not contra- 
dict the formulation of the second law about the arrow 
of time pointing all the time from the past to the future. 



X. FROM THE SIMPLE MODEL TO 
REAL- WORLD PHYSICS 

I believe that the model of generalized baker map is 
interesting even in itself. However, it was studied with 
the hope that its analysis will shed light on the prop- 
erties of more realistic physical model systems. By the 
latter I mean mathematical models, the essential aspects 
of which mimic the properties of corresponding real phys- 
ical objects. Let us start, therefore, the concluding dis- 
cussion by reconsideration of the basic assumptions used 
in the study of the generalized baker map. 

All my numerical simulations were restricted to ratio- 
nal coordinates. The reason was to circumvent artificial 
irreversibility, which would otherwise appear as a numeri- 
cal artifact due to rounding errors. Rational numbers are 
a dense subset in the field of real numbers, and almost 
everything mentioned above for rational coordinates of 
points will apply to real ones as well. There is actually 
only one major problem which the exclusive use of ra- 
tional numbers could cause in our study: the absence of 
periodic orbits for trajectories containing points with ir- 
rational coordinates. This to me does not seem to be an 
essential problem, if we are aware of it. 

Many "realistic" (in the above sense) Af-particle model 
systems (e. g., dilute gas models with Lennard- Jones [2l| 
or hard-sphere |2_2| potential) share a sensitive depen- 
dence on initial conditions with GBM. B 2 is known to 
be strongly mixing [3] and the same is confirmed for 
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many classical systems : 23J. The mixing property is a 
consequence of Lya pun ov instability, discovered in many 
interesting models 22] . The evolution of systems with 
this property leads to local contractivity of phase vol- 
ume (measure) in the forward direction of time (imply- 
ing the opposite when tracking the evolution backwards), 
or to uniform coverage of phase space by the trajectory. 
Approach to an attractor in GBM is then analogous, in 
these more realistic systems, to approach to nonequilib- 
rium steady state, or to equilibrium. 

The notion of absolute age enables us to view the past 
and the future as related to (strange or homogeneous) 
repellors and attractors — the structures in phase space, 
defining the arrow of time. The age, as introduced here, 
relies on the property of B w , that the two coordinates of 
phase points behave differently under B w . It is evident 
that this results from different contractivity of the map 
in different directions. We can therefore expect some 
signs of the absolute age in all systems with strange or 
homogeneous attractors. 

We pointed to the different role of positions and mo- 
menta in the transition to macroscopic description. The 
momenta are responsible for the microscopic difference 
between past and future steady states, but they do not 
help us in their macroscopic identification. This state 
of affairs lies at the basis of macroscopic indiscernibil- 
ity between past and future and thus supports our view 
that macroscopic description is intimately connected to 
the loss of information about the system under consider- 
ation This may best be seen in the fact that we are 
readily able to impose macroscopic constraints on posi- 
tions, but we are not able to do the same (or at least not 
to a comparable extent) with momenta. This dramati- 
cally reduces our ability to influence the states of physical 
systems, thus leading to apparent irreversibility. 

In this sense, the GBM manifests features which co- 
incide with the general properties of e. g., rarefied real 
gases, as regards the global behavior in the distant past 
and future — as far as the interrelation between re- 
versibility and irreversibility is concerned. The "mecha- 
nisms" by which irreversibility appears in reversible sys- 
tems are common to GBM and strongly mixing mod- 
els of TV-particle systems. We may thus say that GBM 
represents an appropriate — albeit extremely simple — 
model, demonstrating essential aspects of interrelations 
between reversibility and irreversibility. The properties 
which became apparent in the study of GBM are there- 
fore applicable to more realistic models, as follows. 

In real experiments, we prepare a state by some manip- 
ulations with particles. The mathematical counterpart of 
such preparation is the localization of phase point ujq of 
a realistic system in its phase space. The evolution of the 
state then results from physical laws and is represented 
as the action of mathematical operator St. Whatever the 
physical preparation of a state may be, the localization 
will always result in a set of numbers having absolute pre- 
cision. Of course, we do not know any state with infinite 
precision, but from the viewpoint of classical mechanics, 



any state must be absolutely precise per se. Application 
to loq of St, simulating relevant physical law, yields ab- 
solutely precise numerical results for u> t , for any value of 
t > 0, provided the system is deterministic. 

As far as GBM is concerned, one can observe (on a 
monitor) even evolution going from repellor — i. e., from 
the distant past — to the present. Such evolution will 
take only finite number of iteration steps: after reaching 
wo, the point will move (forever) to the attractor. How- 
ever, we cannot reconstruct past evolution in real phys- 
ical experiments, since we cannot prepare a state corre- 
sponding to the precise localization of the point Tuj t . A 
notable exception is provided by experimental manipu- 
lations such as those used in spin-echo experiments |24| . 
mimicked in numerical experiments by mathematical re- 
versal of future states. However, due to the sensitive 
dependence of uj t on luq, the slightest deviation from the 
exact value of Tu>t would cause a completely different 
trajectory. 

The same applies to the preparation of states accord- 
ing to a distribution function. In both cases, we pre- 
pare states and functions which are "present" and we are 
not able to prepare "past" states nor functions, which 
would eventually evolve into the present ones. Our pos- 
sibilities to prepare states of physical systems are rather 
crudely restricted to states which are in the absolute 
sense "present". 

In my view, this explains the seeming contradiction 
between microscopic reversibility and macroscopic irre- 
versibility as rooted already at the microscopic level. Or 
still in other words, the problem of irreversibility has pri- 
marily nothing to do with the micro/macro dichotomy. It 
is connected to the past / future dichotomy instead, man- 
ifesting itself in the asymmetry of preparation of states 
close to the repellor/attractor. The macroscopic level 
just adds a further restriction, consisting in our impos- 
sibility to observe or practically influence momenta of 
particles. 

The view presented in the paper may seem to contra- 
dict the orthodox approach, stressing the role of a large 
number of particles in obtaining macroscopic physical be- 
havior, but this role is, actually, not disputed here, it 
is just not given the utmost significance. One certainly 
recognizes that the possibility to prepare a "past" state 
depends crucially also on how large the system is. More- 
over, what is claimed, e. g., in [2£j, is that "the equality 
between the physical entropy production and the dynam- 
ical phase space contraction rates of Hamiltonian ther- 
mostatted systems only holds for macroscopic systems" . 
Restricting the scope of the paper to the investigation of 
irreversibility only, such arguments do not apply to it. 



XI. CONCLUSIONS 

We have analyzed the parametrized set of generalized 
baker maps as simple models of time-reversible evolu- 
tion, exhibiting signs of irreversibility. The analysis has 



15 



shown that it is possible to define past, present, and fu- 
ture states and distribution functions on the definition 
domain of the map. The question about the origin of the 
preferred orientation of evolution the arrow of time 
- can then be answered by stating that it is rooted in 
contractivity of the map: the latter produces a unique 
direction of evolution which can be brought into relation 
with the arrow. This applies equally to the measure- 
preserving classical baker map. 

In an effort to view the map as a model — albeit a very 
crude one — of real systems, and to apply the results of 
the model to them, one sees that a kind of correspondence 
can be established between GBM and macroscopic sys- 
tems, the dynamics of which has an attr actor /repcllor 
pair (visible or hidden). The second law of thermody- 
namics for such systems can then be explained in the 
same way as the unique direction of evolution in GBM. 
Namely, it can be regarded as a consequence of the dy- 
namics which — if continued to encompass the past - 
would lead from the past repcllor to the future attrac- 
tor. The observed properties of the systems then result 
from this simple symmetry. Macroscopic observation of 
systems — the sole viewpoint available to us — hides 
this fundamental microscopic property and causes a se- 
ries of seeming contradictions between the two directions 



of time. Accepting this standpoint, we can claim that 
there is no real "problem" with the law and that the ob- 
jections posed against it (including the classical ones of 
Loschmidt and Zermelo) are simply the results of misun- 
derstanding. 

I think that the described model essentially explains 
the origin of observed (macroscopic) irreversibility in 
unobservably (microscopically) time-reversible systems. 
The observed macroscopic irreversibility is actually re- 
versible in principle, or the microscopic reversibility is 
observationally irreversible. The driving force of such be- 
havior is a sensitive dependence of the evolution on initial 
conditions, which represents the hallmark of chaos. 
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